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Abstract 

We have observed specular reflection of positro- 
nium, Ps, and established that there is ad- 
equate intensity at higher energies to make 
further study worthwhile. The scattering ap- 
pears to be restricted to the outermost surface 
with a mean free path of (0.75 ± 0.15)A for Ps 
in LiF(lOO). With a greater intensity Ps beam 
one should see higher order diffraction beams 
as the result of the periodicity of the surface, 
Ps diffraction thus offers the possibility of be- 
ing a novel and valuable probe to study the 
outermost surface and to study adsorbants on 
it. Two methods for producing Ps beams are 
described. 


Introduction 

At the Brookhaven National Laboratory we 
have initiated a program to study the interac- 
tion of positronium, Ps, with the surfaces of 
solids. An experimental investigation of the re- 
flection of positronium (Ps) from solid surfaces 
is well warranted because of the fundamental 
nature of the electronic processes involved in 
Ps reflection and the possibility of developing 
a valuable new tool for surface structure de- 
termination. Our program of investigating Ps 
reflection was initiated as a result of the fol- 
lowing reasoning. Since Ps can normally be 
expected to undergo elastic collisions from only 
the outer atomic layer of a solid, low energy Ps 
diffraction (LEPSD) could be a unique probe 
of ordered surface structures. This is some- 
what similar to the situation for helium atom 
diffraction, 1 which is a powerful tool in surface 


structure determination because it is only sen- 
sitive to the outer surface layer. However, the 
savings in complexity of He atom diffraction 
by not having to treat multiple scattering from 
subsurface layers, as in the case of low energy 
electron diffraction 2 "" 4 (LEElJ)is somewhat mit- 
igated by having to deal with long range forces 
that dominate in the diffraction. The ~ O.C^eV" 
energies necessary for He atoms to have ~ lA 
de Broglie wavelength results in the He atoms 
having classical turning radii far enough from 
the individual ion cores that the main scatter- 
ing is due to the average potential presented by 
the surface. 1 This requires an accurate treat- 
ment of the atom-surface interaction potential 
which is difficult to obtain and is further com- 
plicated because the depth of the van der Waals 
attractive potential at the surface is approxi- 
mately the same as the kinetic energy of the 
helium atom. 5 In order for Ps to have « lA de 
Broglie wavelength, its energy must be on the 
order of « 75eV\ At this energy, Ps "atoms” 
would be oblivious to the mean surface poten- 
tial and only undergo elastic reflection in close 
encounters with the ion cores. Because of the 
large break-up probability of Ps (binding en- 
ergy = 6.8eV), multiple scattering and other 
subsurface contributions to the elastically scat- 
tered Ps are expected to be negligible. Thus, 
Ps diffraction offers the possibility of being a 
valuable probe. 6 

The degree to which Ps scatters only from 
the outer surface layer is determined mainly 
by the interstitial density of valence or conduc- 
tion electrons of the material. Because of the 
low mass of weakly-bound electrons, and hence 
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large recoil, collisions with the electrons destroy 
the coherence of the scattered Ps and thus must 
be regarded as a source of attenuation of the 
incident and diffracted Ps beam. Typical elas- 
tic cross sections in the lOeV region for Ps-free 
e~ collisions are on the order of 10“ 15 cm. 7 ’ 8 
Thus for solids having interstitial electron den- 
sities of ^ 10 23 cm~ 3 , 9-12 a mean free path 
of ^ lA for the Ps can be expected. Con- 
sequently, LEPSD from a solid surface would 
yield diffracted Ps intensities versus incident 
energy (i.e., ’T(V)” curves) which would be 
dominated by the elastic scattering from only 
the outer layer atomic distribution. In the case 
of an ordered adsorbate overlayer, chemisorbed 
to a surface, however, the incident Ps could eas- 
ily penetrate the relatively open spaces between 
the adsorbate atoms. This would lead to the 
interesting case of interference between Ps scat- 
tering from the adsorbate and from the outer 
surface with a high sensitivity to the structure 
of the adsorbate layer and outer surface. 

In this paper we will first discuss the re- 
sults we have already obtained, 13 present an 
interpretation of them, then provide a descrip- 
tion of the method (gas cell) used to produce 
the Ps beam, and conclude with a description 
of a completely different method to obtain a Ps 
beam which is presently under construction. 

Positronium Specular Reflection from LiF 
Upon entering the experimental chamber the 
Ps beam divergence is limited by aperture to 
5° full width of half maximum of the peak. It 
is reflected from the sample, S, shown in Fig. 1 
and detected by its annihilation gamma rays 
by two bismuth germanate (Bi 3 Ge 40 12 , titled 
BGO) detectors in coincidence. The incidence 
and reflected angles Oi and 0 r are measured 
with respect to the normal to the sample; the 
total angle with respect to the Ps beam is 
~ 0i + 0 r - The sample, S, can be rotated 
through an angle of 50° < 0{ < 90°. The 
intensity of the detected Ps beam from the 
gas cell is measured by removing the sam- 
ple and placing the annihilation plate and 


BGO detectors at ^ — 180°. The distance 



Figure 1 — The experimental arrangement: gc-gas 
cell, s-sample, t-tubes and grids, and a-annihilation 
plate. 

from the center of the gas cell to the annihi- 
lation plate and BGO detectors is a constant, 
40cm, for 100° < ^ < 180°. The efficiency for 
producing and detecting Ps is shown in Fig. 2. 



Figure 2 — The efficiency of detecting Ps at the 
annihilation plate versus the Ps energy. This efficiency 
was obtained with the gas cell filled with Ar at a 
pressure of 10 -s iorr. The Ps atoms were restricted to 
a 5° cone (FWHM). The efficiency includes the decay 
in flight and the Ps detection efficiency. 
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The efficiency is measured per incident positron 
versus the Ps energy. It reflects the efficiency 
of the gas cell filled with 10 ~*torr Ar to form 
ortho Ps in a 5° cone, the decay in flight 
of the initial Ps beam, and the ratio of the 
efficiency to detect of the annihilation gamma 
rays from Ps to those from positrons. The 
inset in Fig. 2 is an expanded scale of the 
low energy Ps region. The absolute efficiency 
is not required for the measurement of the 
reflection coefficient because the efficiencies of 
Ps formation, reflection, and detection are the 
same (except for 2S Ps) for detecting Ps with 
the sample removed from the beamline ( ip = 
180°). We calculate the ratio of the reflected 
Ps intensity to the intensity at ij) = 180°. 
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Figure 3 — Reflection probability of Ps at constant 
incident energy and toted scattering angle jp — 0» + &r 
vs the incident angle, 0*. 


Rocking curves were obtained of the spec- 
ular reflection of Ps from LiF(lOO) by rotating 
the LiF crystal with respect to the Ps beam 
and holding the detectors fixed. Data was 
taken with the position of the Ps detector at 
if; = 100°, 120°, and 130°. The results for 
three different Ps energies and three different 
specular angles are shown in Fig. 3. 13 It is ev- 
ident the Ps specular reflection does occur and 
that the reflected fraction is surprisingly high 
(30 ± 5)% at a Ps energy of 7eV. We also 
measured the fraction of Ps reflected at a fixed 
specular angle as a function of the energy of the 
Ps (see Fig. 4). The fraction of Ps reflected for 
i/> = 100 ±4° from LiF(lOO) was measured when 
the LiF was at a temperature of 160 ± 10°C 
and at a temperature of 300 ± 10°C\ The two 
measurements are in good agreement, although 
the reflection at 300 ° C is somewhat higher than 
that at 160°C. 



Figure 4 — Ps reflection coefficient as a function of 
the Ps energy for = B r = 50.4°. 

A simple interpretation of the Ps reflec- 
tion fraction and its energy dependence can be 
obtained by considering plane waves reflecting 
from a potential step. Letting z be the coor- 
dinate perpendicular to the crystal surface we 
consider the potential: 

V(z) = 0 for z < 0, i.e., in the vacuum out- 
side of the crystal; V(z) — V r + iVim f° r z > 


139 





i.e. inside the crystal there is a real and imag- 
inary potential. Using the one dimensional 
Schrodinger equation in its time independent 
form we obtain the reflection probability: 


R = R 0 


kfj I'l 

k a + ki 


( 1 ) 


where and k Q are the perpendicular compo- 
nents of the Ps wave vectors inside and outside 
the crystal, and where the factor R a is inserted 
to account for the reflection probability is less 
than unity at low energies. (This may be due 
to only a fraction of the surface being clean 
enough). We estimate V r = 4eV, which is the 
difference between the binding energy of Ps in 
the vacuum state (« 7eV ) and the binding en- 
ergy of Ps inside LiF (3eF). 14 If we consider 
= 0 and V r — 4eV we do not predict the 
observed energy dependence of the reflection 
probability. R equals R Q for energies less than 
V r /co$ 2 0i j but it falls off sharply for higher 
energies, and finally approaches approximately 
a 1/Ep s dependence (see the dashed curve in 
Fig. 5). Choosing a different V r only causes a 



Figure 5 — Ps reflection probability vs the incident 
Ps energy. The solid line is calculated with the use of 
Eq. (2) with V r = 4eV and A taken from the fitted line 
in Fig. 6. The dashed line is calculated with V r = 4eV 
and A = oo. 


the data. There is evidently much more appar- 
ent elastic scattering than can be accounted for 
by the real part of the inner potential alone. We 
can obtain a better fit to the data by adding 
an energy-dependent imaginary part, V im , to 
the potential. The wave vector inside the crys- 
tal has two components fc, n = k r + ikim * We 
can solve Eq. (1) for the value of Ar,* m given a 
certain reflectivity R at a given energy E : 


2 2m E 

c ={Ro + R)/(Ro - R) 



where E — Ep 3 cos 2 0i . By using the data 
in Fig. 5, with R = 0.30 dt 0.05 and letting 



E*- E.* E |h feV) 


Figure 6 — Ps mean free path calculated from Eq. (2). 
The straight line is a least-square fit of A = A 0 -f aEp, 
to the data in the interval 16.5eV < Ep s < 56.7eV\ 
The fitted parameters are A — (0.57 db 0.06)A and 
a ~ (4.4 =h 1.7) x 10~ 3 AeV’~ 1 , with a P er degree of 
freedom x 2 / n ~ 15.26/26. 

V r == 4eV we can calculate The mean free 
path is: 


translation of the dashed curve along the hor- 
izontal axis, it does not provide a better fit to 
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where ^ = Aexp ( ifcj n z ) — Aexp i(k r + i&j m )z. 
The result plotted in Fig. 6 is A = (0.75±0.15)A 
for 16.5ey < Ep s < 56.7eV . Below 16.5eV 
the mean free path by this analysis becomes 
unphysically small. The solid line in Fig. 6 is 
a two parameter fit to the data that suggests 
A is slowly increasing with energy. The solid 
line in Fig. 5 is the corresponding reflectivity 
calculated with the use of Eq. 2; the dashed 
line was calculated with a real potential only, 
i.e. A = oo and V r = 4eV. The reasonably 
high-elastic Ps reflection probability observed 
even at high energies, in spite of the presence 
of the large absorptive potential is in retrospect 
not surprising in view of the requirements of 
unitarity. 15 The reason for the few percent 
reflection at energies much greater than V r is 
principally the short A. If Vr were to vanish 
there would still be a measurable reflection 
coefficient for small A. In addition one should 
consider the region between lOeV and lQ.5eV 
may exhibit a higher reflectivity due to a elastic 
scattering from the outer most ion cores. It 
is interesting to note that since the lattice 
parameter of LiF is 4. 02 A, at an incident angle 
of 50° the first order Bragg diffraction would 
occur at a Ps energy of 21eV; however due to 
the short mean free path at low energies it is 
hard to envision the high reflectivity would be 
due to Bragg diffraction at this energy. 

It is evident from the above mentioned 
measurements and analysis that there is a high 
reflection coefficient and short mean free path 
for Ps in LiF, and that an intense, well col- 
limated, and monoenergetic Ps beam holds 
promise as a unique probe of surfaces. 

Positronium Beam— Gas Cell Production 
Low energy positrons emitted from Curare 
magnetically transported through an E x B fil- 
ter out of the shielding blockhouse 16 into an 
Ar gas chamber (see Fig. 7). The pressure of 
Ar in the gas cell was kept at 10 3 torr. By 
the use of baffles and by differentially pump- 
ing the pressure in the experimental chamber 
was reduced to 10 ^tovT, This relatively high 
pressure may not have had much effect on the 


cleanliness of our sample because prior to in- 
troducing Ar, which had a purity of 99.995%, 
into the gas cell; the pressure in the experimen- 
tal chamber was 10~ 10 torr . Thus the gas in 
the experimental chamber was primarily due to 
Ar from the gas cell. Ps is formed in the gas 
cell by the positron picking up an electron from 
the Ar atom. 17-19 The ionization potential for 
Ar is 15.8^, however, the binding energy of 
Ps is 6.8eV thus the threshold energy for Ps 
production is a positron energy of 9.0eV\ The 
first excited state of the Ar atom is ll.SeV 
and the first excited state of Ps is 5.1eV above 
the ground state. Thus a positron beam which 
an energy between 9.0eV and 14.1 eV will pro- 
duce a monoenergetic beam of Ps in the energy 
range of 0 — 5.1eV. The ratio of the cross sec- 
tion for e - ^ to produce Ps in the 2S state to 
producing Ps in the IS state reaches its max- 
imum value of 13% for an e~^ energy of 50eV 
in He. 20 (We are unaware of any calculations 
for a similar ratio in Ar). Above 20.5eV it is 
energetically possible to produce Ps after ex- 
citing an Ar atom. However, the probability 
of a positron undergoing both collisions in the 
gas cell is exceedingly low. In summary, we es- 
timate that the excited state contamination of 
our beam to average less than 5% and the en- 
ergy purity due to Ar excitation to be less than 
1 %. 

Positronium Beam— Foil Production 
We are presently constructing a new high in- 
tensity positron beam in the Material Science 
Building across the street from the High Flux 
Beam Reactor (HFBR) Building. The new fa- 
cility will have four advantages over the present 
one located at the reactor: the low energy 

positron beam will be more intense than the 
beam used for our first generation Ps exper- 
iment described here; we will not be subject 
to the increasingly more severe security regula- 
tions which exist on the operating level of the 
HFBR; the background radiation will be greatly 
reduced; and the area available for experiments 
will be increased. The new facility will have a 
blockhouse approximately twice the size of the 
one built in the reactor building. It is planned 
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Figure 7 — Schematic view of the apparatus to produce a positron beam and a Ps beam by the gas cell 
method. 
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Figure 8 — Schematic view of the apparatus to produce a positron beam and a Ps beam by the foil method. 







to have three low energy positron beam ports, 
the present blockhouse has one beam port. Two 
ports (ports B and C in Fig. 8) will have a elec- 
trostatic transport system and the third (port 
A) will be a magnetically guided transport sys- 
tem. A copper pellet will be irradiated in the 
core of the HFBR for 2 days, then removed 
from its capsule in a blockhouse located along- 
side of the reactor on the operations level of the 
HFBR and deposited into a lead and heavi-met 
container. It will then be transported by an 
electric truck to the Material Science Building, 
and inserted into a crucible in the new block- 
house. The copper pellet will be evaporated 
onto the inside of a cone which then will have 
deposited on it 10 4 A of solid Ne. Our past 
moderator was crystalline copper, but our fu- 
ture one will be solid neon because of its higher 
efficiency. We have measured an efficiency of 
^ 1% for producing low energy positrons with 
solid neon in a cone configuration. 21-24 



Figure 9 — The energy resolution for observing Ps in 
a pulse of 8n* width vs distance from the foil to the 
detector for various Ps energies. 

Instead of using the charge transfer of an 
electron from the Ar atom to combine with a 
positron to produce Ps we will send positrons 
through a thin carbon foil to produce them. 25 
The positron beam will be pulsed which al- 
lows us to perform time of flight measurements. 


This method has three distinct advantages in 
comparison to the gas cell method to produce 
Ps. The advantages are first the surface of the 
sample will not be as easily subject to contami- 
nation as it is from the gas cell since we will be 
able to operate in a ^ 10“ 10 <orr environment. 
The second advantage is that we will be able to 
have a direct measure of the energy of each Ps 
atom with better resolution. Third, although 
the Ps produced by the foil would not be as 
monoenergetic as can be obtained using a rare 
gas target the time of flight method would give 
us the advantage of being able to investigate 
many energies at once. 



Figure 10 — The energy resolution for observing Ps 
in a pulse of O.ln* width vs distance from the foil to 
the detector for various Ps energies. 

The new Ps beam will be produced by the 
following method. 26 ’ 27 The positron beam will 
be transported out of the blockhouse electro- 
statically. It will then enter a magnetically 
guided section. It is then remoder ated to 
minimize its transverse energy component and 
injected into a pinched magnetic field to enter 
a magnetic bottle. Upon entering the bottle 
the beam transverses a rf cavity which oscil- 
lates at 430MHz to give it transverse motion 
and to excite the positron cyclotron resonance. 
The positrons are then reflected by a positive 
potential on a grid, they again transverse the 
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reso- 


rf cavity and they now have enough transverse 
motion to be reflected by the pinched mag- 
netic field at the entrance of the bottle. The 
positrons now oscillate between the two ends of 
the trap for lOOjis. During this time positrons 
continue to stream into the trap from the ra- 
dioactive source and accumulate in the trap. 
After a 100/z.s, the grid which was at a positive 
potential is given a negative potential pulse 
by an oscillator operating at 10kHz and the 
positrons enter a section of 100 ring electrodes 
each having a bias such that the accumulated 
positrons experience a potential which vanes 
harmonically with distance along the beamhne 
and are bunched to a pulse of width 8ns. The 
total length of the accumulator/ pulser/buncher 
is 90 inches. This accumulator/pulser/buncher 
has been constructed, tested and measured to 
have an efficiency of 63%. 20 Upon leaving the 
buncher the positrons are again remoderated 
and enter a second buncher of length one inch 
which produces a harmonic potential due to 
a geometric distortion of the electric field at 
one end of the buncher. This second buncher 
will reduce the positron pulse width from 8ns 
to a subnanosecond bunched positron beam. 
The technical details of the first accumula- 
tor/pulser/buncher are given in Refs. 26 and 27 
and the details of the second buncher are given 
in Ref. 28. Upon leaving the second buncher, 
the positrons traverse a carbon foil and produce 
a pulsed beam of Ps atoms which then enter 
the experimental chamber which was described 
above. The charged particles in the beam are 
removed by electric fields. The energy of the 
Ps atom is determined by measuring the time 
difference between the time the pulsed positron 
beam strikes the carbon foil and the time the 
positronium atoms travel a fixed distance to 
the Ps detection system. The pulse width of 
the positrons leaving the buncher very much af- 
fects the measured energy resolution of the Ps 
atom. This is shown in Figs. 9 and 10 which 
plot the Ps energy resolution versus Ps flight 
path for various energies of the Ps atom. Fig. 9 
is for a positron pulse width of 8 ns and Fig. 10 
is for a positron pulse width of 0.1ns. Although 


a longer flight path increases the energy 
lution it also reduces the number of Ps atoms 
detected due to Ps decay. The effect of Ps de- 
cay in flight is given in Fig. 11 which plots the 
Ps attenuation coefficient as a function of Ps 
energy for various length of flight paths. 



Figure 11 — The attenuation coefficient for Ps 
decay in flight vs the Ps energy for various flight path 
distances. 

The table below shows the expected effi- 
ciencies of the various components of the pulsed 
Ps beam production. The average efficiency was 
calculated for the production of Ps in the en- 
ergy range of GeV to lOOeU in a solid angle of 
10 -2 steradian after traveling a path length of 
20cm. The column at the right lists the num- 
ber of particles to be expected at each stage 
of the beamline. These numbers are predi- 
cated on a spherical copper pellet of weight of 
0.86g and diameter 0.57cm placed in the core 
of the HFBR for a period of 48 hours where the 
positron-emitting isotope 64 Cu is produced by 
the reaction 03 Cu(n/y) 04 Cu. The activity of 
the pellet after 48 hours is 100.5Ci of positron 
emission . 10 Our past experience indicates that 
approximately one-third of this activity can be 
evaporated unto a surface because of losses due 
to decay during the time period from removing 
the source from the reactor to completing the 
evaporation of the Cu on a surface, and due to 
Cu vapor escaping and not being deposited on 
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the surface. We have obtained a moderating ef- 
ficiency of 1.2% with a solid Ne moderator, 24 
but for the purpose of these calculations we are 
estimating an efficiency of 0.7%. We expect the 
efficiency of the 90° bender in the blockhouse 
to be 90%. The number of low energy positrons 


entering the remoderator given in the table, 
7.8 x 10 9 sec, results from these calculations. It 
is the number of particles during the beginning 
of a run, for 64 Cu decays with a half-life of 12.8 
hours. At the end of a two day run the number 
will be reduced by a factor of ^ 16. 


Ps BEAM EFFICIENCY 


Process 

Efficiency 
in % 

Number of Particles 
e+ / sec 

Slow positrons into 
remoderator 


7.8 x 10 9 

Remoderator 1 

30 

2.3 x 10 9 

Accumulator / pulser /buncher 

63 

1.4 x 10 9 

Remoderator 2 

30 

4.2 x 10® 

Ps formation in foil *) 

8.9 x lO" 3 

3.7 x 10 4 

Loss in three grids f) 

72 

2.7 x 10 4 


*) energy range 6eV to lOOeY in a solid angle of 10 ~ 2 ster. 

after traveling a path length of 20cm. 25 
f) 90% transmission grids 


An extensive research program is planned 
to vary the parameters of the carbon foil to 
improve its efficiency to produce Ps (foil thick- 
ness, positron energy, coatings on foil, other 
foils etc.). Another consideration which will 
be examined is the effect of a 34 Ci of positron 
emission source on the Ne moderator. Will it 
cause the surface to be charged? Will it pro- 
duce a large number of defects in the solid Ffe7 
Both effects could adversely effect its efficiency. 
In the event that Ne moderation is adversely 
affected by a strong radioactive source we plan 
to revert back to a transmission moderator to 
produce low energy positrons and to compen- 
sate for its lower efficiency by using a stronger 
64 Cu source. 

Summary 

Positronium reflection from a surface shows 
promise to be an extremely sensitive surface 
probe. We have shown that the mean free 
path in LiF(100) of Ps atoms in the energy 
range l&.beV < E < 56.7eV is 0.75 ± 0.15A. 13 
This assures us that analysis of Ps diffracted 


intensities will only have to take into account 
the outer most atomic layer, and thus avoids 
the complication encountered in LEED where 
the mean free path is typically an order of 
magnitude larger. 29 

A description is given of two different meth- 
ods for producing a Ps beam. Each method 
has its own set of advantages. The gas cell 
method will produce a monoenergejtic Ps beam 
at low energies, whereas the foil beam will 
allow the sample to be in an ultra high vac- 
uum environment and the energy of each Ps 
atom can be measured by a time of flight tech- 
nique. The work was supported in part by the 
National Science Foundation (Grant No. DMR- 
8620168), and in part by the D ivision of Mate- 
rial Sciences, U.S. Department of Energy, under 
contract No. DE-AC-76CH0016. 
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